Running Coupling in the SU(2) Lattice Gauge Theory 



0. Borisenko 1 , M. Gorenstein 1,2 ' s and A. Kostyuk 1 
1 Bogolyubov Institute for Theoretical Physics, Kiev, Ukraine 
2 School of Physics and Astronomy, Tel Aviv University, Tel Aviv, Israel 
3 Institute for Theoretical Physics, Goethe University, Frankfurt, Germany 



On 
On 



Abstract 



Ch . Scenario according to which the SU(2)-gluodynamics is a theory with a 

—5 | nontrivial fixed point is analyzed from the point of view of the modern Mante- 

ls ' Carlo (MC) lattice data. It is found that an assumption of the first order fixed 

point g = gf of the beta function f3f(g) has no contradictions with existing MC 
lattice data. The beta function parameters are found from the requirement of 
constant values for critical temperature T c /A£ p and string tension ■y/a/Af F 
| in MC lattice calculations at 4/g 2 > 2.30. 

o 

The Monte-Carlo (MC) lattice simulations are one of the main sources of the non- 
perturbative results in the gauge field theories. For SU(N) gluodynamics on lattices of size 
iV T x N 3 MC results are the dimensionless functions of a bare coupling constant (3 = 2N/g 2 . 
The transformation of these functions to physical quantities are done by multiplying them 
with a lattice spacing a in the corresponding powers. The length scale L (V = L 3 is the 
system volume) and the temperature T are given by 



X 

5-S 



N a a, T = (N T a)' 1 . (1) 



The lattices with N a >> N T correspond to the finite temperature models, whereas those 
with N a ~ iV r are identified with zero temperature limit. Requirement N T >> 1 is necessary 
to avoid lattice artifacts. 

The best studied quantities in the MC simulations of pure SU(N) gauge theories are the 
dimensionless string tension (y/aa)MC and the critical coupling (3^ c of the deconfinement 
phase transition. The values of f3^ IC were found for the finite lattices and the extrapolation 
to spatially infinite volume ('thermodynamic limit') N a — > oo has been done (see Ref. |! 
and references therein). In what follows we discuss SU(2) gluodynamics. The MC values of 
the critical couplings f3^f c for different N T and of {\faa)MC for different (3 are presented in 
Table I and Table II, respectively. The data are taken from Ref. ]![]. To define their physical 
values one needs a connection between the lattice spacing a and the bare coupling constant 
g. Such a connection is given in terms of the beta function (3f{g) through the equation: 

m = - 4 a ■ (2) 

The conventional perturbation theory gives the following expansion of the beta function 
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Pf( g ) = -W-^ + otf), i, = iH, 6 1 = |(^) . (3) 

N refers to the group SU(N). Differential equation (fj) with (3f F (g) (§) leads to 

aAf = exp (--^) (b g 2 )- bl/2b « - R(f) (4) 



where Af F is an integration constant of Eq. (^|). Eq. (f|) is known as the asymptotic freedom 
(AF) relation. 

Using Eqs. ([!]) and (|) one can calculate 

and 

The values of T c /Af F (see also Ref. |T|) at different A" T are presented in our Table I and of 
\faj Kjf for different couplings (3 in Table II. One observes a rather strong dependence of 
T c j k^f on N T and y^/Af^ on (3. It means that the perturbative AF relation (|J) does not 
work even on the largest available lattices. This fact is known as an absence of the asymptotic 
scaling. In contrast to the problem with an asymptotic scaling the scaling has been observed 
for the ratios of different physical quantities calculated from the lattice expectation values. 
MC data from Tables I and II give for different lattices almost a constant ratio Q: 

= 0.69 ± 0.02 (7) 

if MC values for T c / Af F @ and sfoj Af F (H) are calculated at equal coupling constants (3 in 
the region (3 > 2.30. It suggests a possibility of the universal asymptotic scaling violation: it 
has been proposed in Ref. that a deviation from the asymptotic scaling can be described 
by a universal 'non-perturbative' (NP) beta function, i.e., (3^ p (g) is the same for all lattice 
observables and it does not depend on the lattice size if N a and N T are not too small. 
The following ansatz was suggested 0: 

aAf = X(g 2 )R(g 2 ) , (8) 

where R(g 2 ) is given by Eq. (||) and \(g 2 ) is thought to describe a deviation from the 
perturbative behaviour. The equation (|) has been expected at g — > so that an addi- 
tional constraint, A(0) = 1, has been assumed. The values of T c /A^ p and a/o"/A^ p can be 
calculated then as 

TcMr = N T X(g 2 )R(g 2 ) ' (9) 

^ (10) 
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A simple formula for the function \(g 2 ) was suggested @: 



AG? 2 ) = «'x P ^] . Ill) 



Parameter C3 in Eq. (0) and a new one, T*/A^ p = const, were considered as free parameters 
and determined from fitting the MC values of T c /A^ p ([]) at different N T to the constant 
value T*/A^ p . This procedure gives: 

T */\NP = 21.45(14) , c 3 = 5.529(63) • 1(T 4 . (12) 

The numerical values of T c /A^ p ([)]) are presented in our Table I. In comparison to T c / Aj^ F 
one observes much weaker N T dependence of T c /A^ p (Q). They become now close to the 
constant value T*/A^ p (|T2"|). Due to Eq. (|7|) the constancy of T c /A^ p ([]) guarantees an 
approximate constancy of the physical string tension y/a/A^ p fllCf ) with the average value 
\fa* / A^ = 31.56 in the region of coupling constant (3 = 2.3 -j- 2.8 (see Table II). 

In spite of the evident phenomenological success of the above procedure of Ref. [§] the 
crucial question regarding the validity of the perturbative AF relation @ at g — *■ is not 
solved and remains just a postulate. Do the existing MC data rule out any other possibility? 
To answer this question we reanalyze the same MC data using the same strategy as in Ref. @ . 
A principal difference of our analysis is that we do not assume the AF relation (|]) between 
g and a at g — > 0. Instead of this standard approach we check a quite different scenario with 
the fixed point (FP) gf of the beta function. Note that the so-called FP field theory models 
were considered a long time ago [[|. This theoretical possibility has been also discussed in 
Ref. ||. It was demonstrated in Ref. that the precise data on deep inelastic scattering 
do not eliminate the FP model and other tests would be necessary to distinguish between 
AF and FP QCD. 

Let us assume that the beta function of the SU(2)-gluodynamics has a zero of the first 
order at some FP gf. We thus have in the vicinity of this point 

Pf P (9) = -b{g- g f ) . (13) 

From Eqs. (0) and ( pT3[ ) we find then 

a A pp = (g - g f )V , (14) 

with A pp being an arbitrary integration constant of differential equation (||). It follows then 
for the critical temperature 

TJA PP = — — (15) 

and for the string tension 



^ /Ai = ( S - jw> ■ (16) 

Our requirement similar to that of Ref. is to fit T c /A pp (fl"5|) for different N T to a constant: 
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T*/Al y = C T = const , (17) 

where numerical value of Ct is a priori unknown. This requirement leads to the following 
expression for the critical coupling 

* = (C ^>-° + «/ • A = [( C T N T )-* + gif ■ (18) 

In addition we require the constancy of \faj Kjf CP~6f) for different (3. It introduces another 
unknown constant 

v^*/A£ P = C a = const (19) 
and leads to the model equation for y/aa as a function of (3: 

v^a = C CT (s - g f fl h =C a (2/^ - g f f/ b . (20) 

Our fitting procedure for finding beta function parameters gj and b as well as the con- 
stants Ct (0) and C a (]T9|) is to minimize x 2 defined as 
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where (3 C and v^aa are given by model equations ( 18|) and (|20|) , Ap^ IC and A( v / aa) A f C < 
stand for uncertainties of the corresponding MC data. We assume that deviations from the 
formula ([14]) are negligible for (3 > 2.30 and use available MC data from Tables I and II 
satisfying this criterion: in the first sum of Eq. (^1|) we include MC points (3^ IC from Table 
I for N T = 5, 6, 8, 16 and the second sum is carried over $ = 2.3, 2.4, 2.5, 2.6, 2.7, 2.85 
with corresponding MC values of (v^gOmc from Table II. 
The x 2 reaches its minimum x m in = 2-15 at 

g f = 0.563 , b = 0.111 , C T = 3.15 , C a = 4.59 . (22) 

For the set of parameters (|22"D the values of T c /Af p flT5|) and \faj 'Af p are presented in 
Tables I and II. They are constant within the errors^ for different lattice sizes and different 
(3 values for (3 > 2.30. In Figs. 1 and 2 our fits for (3 C fll8D and (s/cra) ( p0|) (with the same 
parameter set ( E2p ) are compared with MC values of (3^f c and (^Imc from Tables I and 
II. 

The standard criterion x 2 < Xmin + 1 HI defines an confidence region of the model 
parameters which gives a rather large variance of gf, b, Ct, C a , but with a very strong 
correlation between them. In Fig. 3 we show the projection of this confidence region to 
the (g f ,b) plane. At the end point A, (g f = 0.288,6 = 0.800), we find C T = 0.174 and 
C u = 0.253. Another end point B, (g f = 0.717,6 = 0.142), corresponds to C T = 8.90 and 



1 The errors are induced solely by the uncertainties of the MC data for j3^f c and (\fcra)MC U1 
each point, and do not include the variance of the fit parameters. 
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C a = 12.98. Physical observables stay almost unchanged under large variance of the model 
parameters due to their strong correlations. We find 



T c /y/a = C T /C a = 0.687 ± 0.005 , (23) 

where the error in Eq. (^) caused by the parameter variations is calculated from the variance 
matrix of the parameters 0. 

From Eq. @ one can easily reconstruct the NP beta function j3^ p {g) 



Pf b - b l9 * + 3c 3 g* 



iNP 



In Fig. 4 we compare different beta functions discussed in our paper: [3f F (g) (f|), Pf P {g) 
([2D and Pf p (g) flT3|). It is amazing that (3? p (g) with NP corrections (§) of Ref. § is 
very close to our straight line /3 FP (g) fll~3] ) in a rather wide region of the coupling constant. 
Note also that different numerical values of T*/A^ p and T*/A FP are caused by a freedom in 
choosing the form of expression ([14]) defining the arbitrary integration constant of differential 
equation (f2|). 

At the present moment there are no SU(2) MC data for the critical temperature and the 
string tension in the region f3 > 2.85. In the framework of our FP scenario we can predict 
the values of \faa and (3 C for the future MC calculations. They are presented in Tables III 
and IV. As is seen from Figs. 1 and 2 these numbers can be hardly distinguished from those 
obtained with NP corrections (|5p @ to the AF relation. One needs very large values of (5 
to observe the difference. The principal difference is of course exist: f3 c — > oo for N T — ► oo 
in the approach of Ref. and (3 C — > 4:/gj = const for N T — > oo in the FP scenario. 

It is clear that obtaining MC data for the critical temperature and string tension in the 
SU (2) gauge theory for the region of j3 essentially greater than 2.85 would require very large 
lattices and hardly possible in the nearest future. It seems that these restrictions are not so 
severe for the finite volume observables in the SU(2) gauge theory. In Ref. the quantity 
g 2 {L) defined as the response of the system in hypercube L x L x L x L to a constant 
color-electric background field was studied for f3 = 2.6 -j- 3.7. To check the consistency of our 
FP scenario with these results we reanalyzed some MC data obtained in Ref. 0. Namely 
we use in our analysis the sets of the bare couplings (3 tuned to achieve constant values of 
fiLj), (j = 0, 2, 4, 6, 8) [] for different lattice sizes N a = N T = N. 

From Eq. fli~4]) it follows that the dependence of g on N at constant g 2 (Lj) is given by 




9j {N) = V +9f, (25) 



where relation a = L/N has been used. In contrast to our previous consideration, the cutoff 
dependence (lattice artifact) is not negligible in the present case. Following to Ref. we 
assume that it is proportional to 1/N. Therefore to fit the data from the Table V we use 
the following formula 



2 We enumerate Lj with successive even numbers to retain the notations of Ref. |fj]. 
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■ 3 ' {N) = (OT - % (26) 



where Cj is a constant. 



In our fit procedure the FP beta function parameters gf and b have been fixed at their 
values ( p2]) found from the previous analysis, and LjA FP and Cj for each column of Table 
V are chosen to minimize the following expression 



The results of the fit are shown in Table VI and in Fig. 5. 

For a comparison we have made the similar fit assuming AF behavior of gj(N), i.e. 
instead of Eq. fl25|) the value of gj{N) was determined as a solution of the transcendental 
equation 



LA AF 
N 



R(g](N)) (28) 



with function R defined by Eq. (|j). As is seen from Table VI, yj/dof is almost two times 
smaller in the FP scenario then that in AF one. This value correspond to the region 1.17 < 
g < 1.25, where the difference between AF and FP beta-functions becomes large (See Fig. 4). 
For the rest of the data sets both approaches give nearly the same fit quality. Still, the FP 
scenario provides a better agreement of the ratios Lj/L 8 with those of Ref. M shown in the 
last column of Table VI. Therefore, MC data of Ref. R are consistent with FP beta function 



behavior (|13|,^). In the region of g ~ 1 we do not observe the difference between AF and 
FP scenarios: as seen from Fig. 4 (3f F and (3j p are close to each other in this region of g. 
Hopefully additional lattice calculations in the spirit of Ref. at bare coupling g < 1 would 
allow one to find the true zero of the beta function (3f{g). 

We conclude that available MC lattice data in the SU(2)-gluodynamics do not exclude 
the possibility of the FP scenario (|T4T). New lattice data are necessary to prove (or disprove) 
the AF relation (IT). 
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TABLES 



TABLE I. MC data for critical couplings /3* at different iV T are taken from Ref. [1]. The 
values of T c j ' Ajf are calculated from Eq. (||) assumed the perturbative AF relation (||) (see also 
Ref. [1]). T c /A^ p are obtained from Eq. (0). Our results for T C /A FP followed from Eq. ( |i~5| ) in the 



FP 


scenario are presented in 


the last column. 






N T 






T C /A N L P 


T c /A pp 


2 


1.8800(30) 


29.7 


8.65(12) 


1.349(16) 


3 


2.1768(30) 


41.4 


18.69(21) 


2.696(29) 


4 


2.2988(6) 


42.1 


21.44(05) 


3.084 (7) 


5 


2.3726(45) 


40.6 


21.95(33) 


3.167(48) 


6 


2.4265(30) 


38.7 


21.81(22) 


3.156(32) 


8 


2.5115(40) 


36.0 


21.44(27) 


3.124(41) 


16 


2.7395(100) 


32.0 


21.50(64) 


3.200(99) 



TABLE II. MC data of {y/aa)MC f° r different lattices and coupling constants are taken from 
Ref. [1]. The values yfa I Ajf and y/a/A^ p are calculated from Eqs.(j6|) and (|i~0l), respectively. 
Last column corresponds to our results for yfa / ' A FP (16) in the FP scenario. 





N T 











v^/Ar 


8 


10 


2.20 


0.4690(100) 


61.7(14) 


28.56 (61) 


4.116(88) 


10 
16 

32 
20 
32 
48 


10 
16 

32 
20 
32 
56 


2.30 
2.40 
2.50 
2.60 
2.70 
2.85 


0.3690(30) 

0.2660(20) 

0.1905(8) 

0.1360(40) 

0.1015(10) 

0.0630(30) 


62.4(5) 
57.8(4) 
53.3(2) 
49.0(14) 
47.1(5) 
42.8(21) 


31.78 (26) 
31.94 (25) 
31.51 (14) 
30.70 (91) 
31.03 (31) 
30.00(143) 


4.574(38) 
4.615(35) 
4.587(20) 
4.509(133) 
4.601(46) 
4.511(215) 



TABLE III. Predictions for (y/aa)MC a t different j3 according to Eq. (p0[). 



p 


2.90 


2.95 


3.00 


3.05 


3.10 


y/aa 


0.0552(14) 


0.0476(15) 


0.0411(16) 


0.03555(17) 


0.0308(17) 



TABLE IV. Predictions for at different N T according to Eq. (|l|). 





20 24 


28 


32 


Pc 


2.8077(49) 


2.8683(71) 


2.9201(93) 


2.9653(115) 



s 



TABLE V. The bare coupling j3 at different lattice size N for fixed g 2 (L) [7]. The uncertain- 
ties of j3 in columns 2-5 were recalculated from errors of g 2 (L), given in Ref. [7], using linear 
interpolation for the dependence g 2 (L) on (3 at fixed N. 



AT 


P 




g 2 (L ) = 2.037 


f{L 2 ) = 2.380 


2 2 (L 4 ) = 2.840 


f(L e ) = 3.550 


g 2 (L 8 ) =4.765 


5 


3.4564(25) 


3.1898(22) 


2.9568(19) 


2.7124(34) 




6 


3.5408(40) 


3.2751(35) 


3.0379(31) 


2.7938(48) 


2.5752(28) 


7 


3.6045(42) 


3.3428(36) 


3.0961(33) 


2.8598(50) 


2.6376(20) 


8 


3.6566(47) 


3.4009(42) 


3.1564(38) 


2.9115(55) 


2.6957(21) 


10 


3.7425(59) 


3.5000(61) 


3.2433(53) 


3.0071(77) 


2.7824(22) 


12 










2.8485(32) 


14 










2.9102(62) 



TABLE VI. The results of fitting the data from Table V with Eq. (26) assuming FP (25) and 
AF (|28| ) scenario. The ratios Lj/Lg are compared with those calculated in Ref. [7]. The number 
of degrees of freedom (dof), the difference between the number of points and the number of fit 
parameters, equals 3 for j = 0, 2, 4, 6 and 4 for j = 8. 





FP 


AF 


Ref. [ 


1 


j 


Xf/dof 


Cj 


Ljk F L p x 10 


Lj/L 8 


Xj 2 /dof 


Cj 


LjA£ F x 10 2 


Lj/L 8 


Lj/Ls 





1.25 


0.315 


0.102 


0.072 


1.08 


0.215 


0.139 


0.080 


0.070(8) 


2 


0.09 


0.559 


0.187 


0.131 


0.09 


0.378 


0.254 


0.146 


0.124(13) 


4 


0.75 


0.422 


0.387 


0.272 


0.81 


0.177 


0.512 


0.293 


0.249(19) 


6 


0.14 


0.585 


0.734 


0.516 


0.18 


0.220 


0.933 


0.534 


0.500(23) 


8 


0.89 


0.589 


1.442 


1.000 


1.67 


0.032 


1.742 


1.000 


1.000 
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FIGURES 




FIG. 1. Circles with errorbars represent the MC data listed in Table I. The bold solid line 
corresponds to Eq. (18) of the FP scenario. The dotted and thin solid lines show f3 c = 4/<£? found 
by solving the equations T*/A^ p = [iV T A(#?).R( ff 2)]--i and T* / kf p = [N T R(g 2 c )}- 1 , respectively, 
with T*/A% p = 21.45. 
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0.250000 - 




0.031250 



0.015625 



FIG. 2. Circles with errorbars represent the MC data listed in Table II. The bold solid line 
corresponds to Eq. (20) of the FP scenario. The dotted and thin solid lines are calculated as 
\(g 2 )R(g 2 )V^/A% p and R(g 2 )V^/A^ p , respectively, with V^/A^ p = 31.56. 
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1.2 1.4 



FIG. 4. 



The behavior of fij p {g) 



33) (bold solid line). 



(5) P 



(dotted line) and (3f F (g) 



(|3]) (thin solid line) . Two vertical lines show the region of the coupling constant (/? 
MC data used in the fitting procedures. 



2.3 



/ 

2.8) of 
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FIG. 5. The results of fitting the data from Table V with the formulas (pq ) and (p6|). 
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